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ABSTRACT 

The re-emergence of the 0.99 ftm FeH feature in brown dwarfs of early- to mid-T spectral type has 
been suggested as evidence for cloud disruption where flux from deep, hot regions below the Fe cloud 
deck can emerge. The same mechanism could account for color changes at the L/T transition and 
photometric variability. We present the first observations of spectroscopic variability of brown dwarfs 
covering the 0.99 /rm FeH feature. We observed the spatially resolved very nearby brown dwarf 
binary WISE JI049I5.57-531906.1 (Luhman 16AB), a late-L and early-T dwarf, with HST/WFC3 
in the G102 grism at 0.8-1.15 /rm. We find significant variability at all wavelengths for both brown 
dwarfs, with peak-to-valley amplitudes of 9.3% for Luhman 16B and 4.5% for Luhman 16A. This 
represents the first unambiguous detection of variability in Luhman 16A. We estimate a rotational 
period between 4.5 and 5.5 h, very similar to Luhman 16B. Variability in both components complicates 
the interpretation of spatially unresolved observations. The probability for finding large amplitude 
variability in any two brown dwarfs is less than 10%. Our finding may suggest that a common but 
yet unknown feature of the binary is important for the occurrence of variability. For both objects, 
the amplitude is nearly constant at all wavelengths except in the deep K I feature below 0.84 ftm. No 
variations are seen across the 0.99 fim FeH feature. The observations lend strong further support to 
cloud height variations rather than holes in the silicate clouds, but cannot fully rule out holes in the 
iron clouds. We re-evaluate the diagnostic potential of the FeH feature as a tracer of cloud patchiness. 
Keywords: binaries: visual — brown dwarfs — stars: atmospheres — stars: individual (WISE 
J104915.57-531906.1, Luhman 16AB) —stars: variables: general 


1. INTRODUCTION 

At the transition from L to T spectral types (T.// ^ 
1200 — 1300 K), the clouds in the atmospheres of sub- 
stellar objects disappear from the photosphere, result¬ 
ing in drastic changes in the spectra over a small 
range of effective temperatures of only 100-200 K (e.g. 
iKirkoatrickl 120051) . The process by which the clouds 
disappear is not yet well understood. A view in 
which increasing partic l e sizes lead to rai n out (e.g. 
iTsuii fc Nakaiima I2003I iKnapp et al.l 120041 ) , and thus 
physically and optically thinner clouds that eventually 
disappear completel y, can broadly explain the color- 
magnitude e volution (iSaumon fc Marlevll2008D and spec- 
tral series (|Cushing et al.l l2008t iStephens et al.l I2009D 
from late L to mid T dwarfs, as well as the luminosi¬ 
ties of an L/T transition b inary with known dynamical 
masses (jPupiiv et al.ll2015l) . 

Another mechanism that has been proposed is th e 
break up of cloud layers (lAckerman fc MarlevI 120011). 
Models of partly cloudy atmospheres (iBurgasser et al.l 
l2?inl IMaHev et a].1IWm can also broadly explain the 
color evolution in the transition, in particular the J 
band brightening. Additional evidence pointing to 
cloud disruption was provided by t he discovery of high- 
amplitude variable early T dwarfs (lArtigau et al.ll200l : 


iRadigan et al.l [201^ iGillon et al.l I2013D . Furthermore, 
the strengthening of the FeH feature at 0.99 ftm for early 
to mi d-T dwarfs has been in terpreted as a sign of cloud 
holes (IBurgasser et al.l[2002i) . Iron condenses into clouds 
and depletes the atmosphere of iron-bearing gases such 
as FeH. The feature weakens for late type L dwarfs, and 
as the clouds sink below the photosphere the detectable 
FeH is expected to decrease further. However, its re- 
emergence for early to mid-T dwarfs suggests that the 
cloud deck does not sink gradually, but opens holes into 
hotter regions in which FeH is not depleted. 

The re-emergence of FeH is based on a fairly small 
sample of brown dwarfs. Recently, observations of 
the L/T transition binary WISE J104915.57-5 3 1906.1 
(hereafter Luhman 16AB) by iFahertv et al.l (l2014f) 
showed that the strength of the FeH feature was 
equal for the L7.5 type component and for the TO.5 
type component, even though the TO.5 component is 
highly variable (IGillon et al.l 120131: iBiller et^ 120131 : 
IBurgasser et al.l[2014D . suggesting patchy clouds. Fur¬ 
thermore, spe ctroscopic variabilit y measurements of 
Luhma n 16B (iB uenzli e t al.l I2015D and other early T 
dwarfs (|Apai et al.l 120131) at 1.1 to 1.7 fim have shown 
that a mixture of thin and thick clouds, rather than 
cloud holes into significantly hotter regions, are required 
to explain the wavelength dependence of the variability 
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Figure 1. A raw first order spectrum on the detector for Luhman 
16A (top) and B (bottom), spanning about 0.8 to 1.15 fim. The 
dispersion direction is approximately toward north east. The gray 
scale is logarithmic. 


amplitude. 

Here we present the first observations of spectroscopic 
variability of brown dwarfs between 0.8 and 1.1 ^m, cov¬ 
ering the 0.99 fim FeH feature. We observed the Luhman 
16AB binary with the Hubble Space Telescope (HST). 
Luhman 16 is an ideal t arget due to its very close proxim¬ 
ity to the Earth (2 pc. lLuhmanll20l^ iBoffin et al.ll20l4 


ISahlmann fc Lazorenkol 120151) . allowing high signal-to- 

noise observations, and the possibility for resolved 
observations of both a late L and early T dwarf. 
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). We detect signifi- 

cant variability in both brown dwarfs. We discuss the 


wavelength dependence of the variability and compare 
to heterogeneous cloud models, with a special focus on 
the FeH feature. 


2. OBSERVATIONS AND DATA REDUCTION 

Observations of Luhman 16A and B were carried out 
during 5 consecutive HST orbits between 17:20 UT on 
2014 November 22 and 00:28 UT on 2014 November 23 
(Program 13640). During each 96 min orbit, the tar¬ 
get was observable for 50 min with gaps of 46 min when 
the target was behind the Earth, except for the first gap 
which was only 30 min. The shorter first gap is a result of 
scheduling near the continuous viewing zone of HST. We 
used WFC3 in its infrared channel, obtaining spectral 
time series of the binary with the G102 grism. The G102 
grism spans wavelengths between about 0.8 and 1.15 ^m. 
We used the 256x256 subarray mode which allowed stor¬ 
age of all observations in an orbit in the WFG3 buffer, op¬ 
timizing the observing efficiency because no buffer dump 
was required during the visibility periods. With a pixel 
size of Ri0'.'13, the field of view was f»30"x 30". The first 
order spectra, spanning ?»155 pixels with a dispersion of 
2.36 to 2.51 nm pixel”^, were fully captured on the sub¬ 
array (Fig. [T]). Several faint spectra of background stars 
are also visible in the field, but none of them overlap with 
any of the spectra of our binary target. The grism disper¬ 
sion direction, given by the orientation of the spacecraft, 
was set very close to perpendicular to the line connect¬ 
ing the binary, minimizing spectral overlap between the 
two objects. For wavelength calibration, a direct image 
in the F098M filter was obtained at the beginning of 
each orbit to measure the location of the sources. All 
observations were conducted in staring mode without 
dithering to avoid errors from pixel-to-pixel sensitivity 
variations. For the spectroscopic time series, we used 
the SPARS25 readout mode, obtaining exposures with 
three non-destructive reads at 0 s, 0.278 s and 22.62 s 
(NSAMP=2). We took 72 exposures in the first orbit 
and 75 in subsequent orbits. The cadence, i.e. the expo¬ 
sure time plus overhead for one exposure, was 41 s. We 


averaged 3 spectra, thus binning to a cadence of 123 s 
to improve the signal-to-noise ratio (SNR). In total, the 
time series spanned 7.1 h. Direct images were taken in 
the SPARSIO readout mode with a read at 0 s and at 
0.278 s (NSAMP=1), two such exposures were obtained 
in each orbit. The maximum number of counts recorded 
in a pixel in the spectra was «20,000, well below the 
regime where image persistence can become significant 
(above ~ 40, 000 counts). 

Data reduction was c arried out in the s a me manner as 
inlBuenzli et al.l (j2012f ) ; lAoai et HI (l2013f ) ; iBuenzli et al.l 
(|20l4 I2015D and readers are referred to these papers 
for details. In short, we used a combination of the 
WFC3 pipeline, custom IDL and Python routines and 
the PyRAF software package aXcQ^ which was specifi¬ 
cally developed to extract and calibrate slitless spectro¬ 
scopic data from HST. The zero-read and dark current 
was subtracted and the images were corrected for non¬ 
linearity and gain, as well as bad pixels and cosmic rays. 
The axeprep routine was used for background subtrac¬ 
tion by scaling a master sky frame, and the axecore 
routine for flat-fielding, wavelength calibration, extrac¬ 
tion of the two-dimensional spectra and flux calibration 
with the G102 sensitivity curve. The extraction width 
for the spectra was set to 6 pixels. We applied aper¬ 
ture c orrection based on the values of iKuntschner et al.l 
(Imh) , interpolating between the listed values. Because 
the beginning of the first orbit is impacted by a well- 
known ramp effect, we remove the first 20 min, as they 
are not crucial to our science resul ts. We note that the 
G141 observations of Luhman 16 in lBuenzli et al.l (I2015D 
had to be taken in SPARSIO rather than SPARS25 read¬ 
out mode because a shorter exposure time was required. 
With the SPARSIO mode we had found a much larger de¬ 
creasing ramp effect spanning the first two orbits. Those 
systematics had been corrected using the non-variability 
of Luhman 16A in those observations. Such a procedure 
was not necessary for the new SPARS25 observations. 

The separation of the binary has further decreased 
within the past 12 months, from 1.24"to approximately 
0.9". We estimate the flux overlap in the extraction aper¬ 
tures of Luhman 16A and B from the aperture correction 
values. In y-direction, the two spectra are separated by 
6.58 ± 0.02 pixels. With a 6 pixel extraction width, the 
flux from A present between 3.58 and 9.58 pixels dis¬ 
tance from the peak on one side is in the aperture of B, 
and vice versa. We find that depending on wavelength 
2.5-3% of the flux of the other brown dwarf is present in 
the aperture, compared to about 90% of the flux of the 
brown dwarf on which the aperture is centered. Consid¬ 
ering that B is about 40% brighter in the Y band peak 
than A and that the maximum variability is less than 
10% for B, less than 0.5% of the measured variability of 
A can be attributed to originate from B. On the other 
hand, only «0.1% of the measured variability of B can 
originate from A. 

3. RESULTS 

We detect significant variability at all wavelengths for 
both sources. The integrated light curves are shown in 
Figure [5J The maximum and minimum spectra and the 
amplitudes as a function of wavelength are shown in Fig- 

^ http://axe-info.stsci.edu/ 
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Figure 2. Light curve for Luhman 16 A (top, red) and Luhman 16 
B (bottom, blue) obtained by integrating counts over the full spec¬ 
tral range of the G102 grism. For compa rison, the J band lig ht 
curve (1.22-1.32 pm) of Luhman 16 B from lBuenzli et al.l Il2016l 'l is 
also plotted (black), shifted in time to match the phase of the dip. 
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ure[3l These observations represent the first unambiguous 
variability detection for Luhman 16A, although tentative 
evidence for vari ability at these wa velengths was previ¬ 
ously reported in iBiller et al.l (|20I3D . 

3.1. Light curves 

Because the variability amplitudes show very little 
wavelength dependence over the whole measured spec¬ 
trum, we only discuss light curves integrated over the 
whole covered spectral region (Fig. [21). The peak-to- 
valley amplitudes are 9.3% for Luhman I6B and 4.5% 
for Luhman 16A. Since we observe a turning point at 
the maximum and minimum flux values, it is very likely 
that the maximum and minimum are indeed covered for 
both objects and we measure the full amplitude. A rapid 
and deep dip in the flux is observed for Luhman I6B in 
the third orbit. This dip is very similar to the one seen in 
the longer wavelength o bservations taken one year ear¬ 
lier (|Buenzli et aLllMTsI) . At other rotational phases, the 
light curve shape has changed significantly. 

With a coverage of only 7 h and significant gaps, as well 
as potential light curve evolution on those time scales, we 
cannot determine reliable periods. We note a potential 
period of Ri5.15 h for Luhman 16B due to the presence 
of two identically shaped local minima in the second and 
fifth orbits. This is slightly lo nger than 4.87 ± 0.01 h 
measured bv iGillon et al.l (120131 ) , but consistent with th e 
5.05 ± 0.1 h perio d measured bviBurgasser et all (j2014D . 
For Luhman IBA. lBiller et al.l (l2013f ) suggest a period of 
3-4 h from their tentative detection of variability. In our 
light curve, we find two maxima, although the second 
maximum is brighter by approximately 2% of the mean 
flux. Assuming a similar light curve shape but overall 
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Figure 3. HST Observations of Luhman 16 A and B with the 
G102 grism. Top: Maximum and minimum spectra averaged over 
min. The location of the FeH and CrH features are marked. 
Bottom: Relative amplitude as a function of wavelength, i.e. the 
difference of the maximum to the minimum divided by the mean. 
The Icr error is shown as a gray lines. 


brightening, we roughly estimate the rotation period to 
be between 4.5 and 5.5 h. Periods below 4.5 h can be 
excluded unless the light curve shape also changes dras¬ 
tically. The period may also be significantly longer than 
7 h in case of a double-peaked light c urve. However, 
theusin? measurement of 17.6 km/s (jCrossfield et all 
iMi) indicates a maximum rotation period of 7 h as¬ 
suming a radius of 1 Rj. Therefore, a double-peaked 
light curve for Luhman 16A is not very likely. The po¬ 
tentially very similar rotation periods for Luhman 16A 
and B, together with different projected equatorial rota¬ 
tional velocities (26.1 ± 0. 2 km/s for B an d 17.6 ±0.1 
km/s for A) measured by (jCrossheld et al.ll2014D . would 
indicate some misalignment between the rotation axes of 
the binary components. 

3.2. Spectroscopic variability 

The wavelength dependence of the variability ampli¬ 
tude is largely flat for both sources. No significant de¬ 
viation in the amplitude is seen in the FeH feature at 
0.99 /im or in most other absorption features. Only at 
wavelengths < 0.83 ^m the amplitude decreases in both 
sources. This coincides with the onset of the very deep 
K I feature centered at 0.77 pm. However, it is also near 
the edge of the spectrum where the throughput as well as 
the flux is low and errors are larger. Additionally, small 
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Figure 4. Principal component analysis of Luhmanl6 A (red) and B (blue). Top Left: First and second PCA components compared 
to the mean flux, normalized to the maximum flux of the B component. The first PCA Components are scaled to allow comparison with 
the mean flux, the second PCA components are scaled for better visibility. Bottom left: Variability amplitude (maximum difference of 
coefficients multiplied by the component and divided by the mean) derived from the mean and first component compared to the variability 
amplitude from the maximum and minimum spectrum (gray, lines from Fig.[3]|. Right: Coefficients of the first and second components. 


dips in the amplitude are seen for Luhman 16B between 
0.85 and 0.9 ^m. 

We estimate the error on the variability amplitude as 
a function of wavelength by calculating the standard de¬ 
viation of the flux in neighboring pixels, which belong 
to one resolution element, after normalizing with a mean 
spectrum over the full time series. While the error is 
below about 0.1% at A > 0.9 fim for both sources, it in¬ 
creases rapidly toward shorter wavelengths (see Fig. [ 31 ). 
At most wavelengths, the SNR of the variability ampli¬ 
tude is between 20 and 50 for Luhman 16 A, and can 
rise above 100 for Luhman 16B. The SNR remains above 
5 down to 0.82 /im for both objects, where the ampli¬ 
tude has decreased by about 1% for Luhman 16A and 
by 3% for Luhman 16B. To test whether the decreasing 
amplitude is not only a result of increasing error due to 
the lower signal, we fit a linear model between 0.835 and 
0.8 /im and compare the chi-squared value to that for a 
model where the amplitude remains constant at a fixed 
value of 4% or 8% respectively. The chi-squared value 
is larger by 25 for Luhman 16A and by 100 for Luhman 
16B. Despite the two additional degrees of freedom for 
the linear model with a slope, we can reject the null hy¬ 
pothesis that the constant model is valid at > 99.999% 
confidence. The decrease in amplitude as a function of 
wavelength therefore appears to be highly statistically 
significant for both objects. This statistical test assumes 
Gaussian errors. We cannot fully exclude a systematic 
effect near the grism edge that would introduce corre¬ 
lated errors. However, no similar amplitude change is 
seen at the other edge of the spectrum. 


The achie ved precision is com parable to that for the 
G141 grism ([Buenzli et al.ll2015h a. except at the shortest 
wavelengths where the flux drops dramatically. This is 
mainly due to the fact that Luhman f6A and B are very 
bright and the G141 observations were less efficient due 
to very short exposure times. For fainter brown dwarfs, 
for which the observing efficiency is comparable in both 
grisms, the Gf02 observations will typically have lower 
SNR due to the overall lower count rate. 

We apply a principal component analysis (PCA) to 
determine the number o f variable componen ts, anal¬ 
ogous to the analysis in iBiienzli et al.l (|201,5l ) for the 
WFC3/G141 observations of Luhman 16B (Fig. [4]). As 
was the case at longer wavelengths, a single variable 
component on top of the mean can explain about 93% 
of the observed variability of Luhman 16B. The higher 
order components can largely be attributed to noise. 
The flux F{X, t) can therefore approximately described 
as F{X,t) Ri {F{X,t))t + ai(t)ei(A), where ai{t) is the 
first coefficient and ei(A) the first principal component 
(eigenvector). 

The first component ei(A) is essentially proportional to 
the mean flux, as expected for a flat amplitude as a func¬ 
tion of wavelength. The reconstructed time series using 
only the first component confirms the drop in amplitude 
at < 0.84 /im and the small dip between 0.85 and 0.9 /im. 
Additionally, a slight gradient at wavelengths < 1 /im is 
revealed. 

The PCA analysis provides similar results for Luh¬ 
man 16A, although the first component accounts only for 
about 66% of the variability. Nevertheless, all other com- 
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ponents account for well below 2% and are therefore also 
mostly noise. Again, the wavelength dependence of the 
amplitude is nearly flat, except for a small dip at about 
0.85 /im and a rapid decrease below 0.83 /im. There is 
a very small gradient in the opposite direction than for 
Luhman 16B, with a ~0.5% decrease in amplitude from 
0.88 to 1.13 ^m. 

4. DISCUSSION 

4.1. High amplitude variability for both components of 

a binary 

The discovery that Luhman 16A is unambiguously 
variable with an amplitude of about 4%, at least some 
of the time, while Luhman 16B is has the second largest 
J band amplitude of any known brown dwarf, poses the 
question whether the Luhman 16 system has any partic¬ 
ular feature that may favor the occurrence of high am¬ 
plitude variability over average field dwarfs. 

We first estimate the statistical probability of finding 
a variable late-L and early-T dwarf assuming indepen- 
den t probabilities. We u se t he combined sample from 
the iRadigan et al.l (j2014[ l and iWilson e t al.| ([20 1^ sur - 
veys, adjusted for the re-analvsis from iRadiganI (j2014[l . 
to derive variability frequencies. Both were large ground 
based photometric surveys in J band. Although our 
variability detection for Luhman 16A is in the Y band, 
the measured wavelength dependence and models sug¬ 
gest that we would expect a similar amplitude in the J 
band (cf. 14.21 and I4.3|l . From all the sources, we pick 
the spectral type range L6-L9 to estimate the variability 
fraction of late-L types like Luhman 16A, and L9.5-T2.5 
for early-T types like Luhman I6B. We remove 5 targets 
with noisy observations where the detection limits might 
not have allowed a confident detection of the variability 
of Luhman 16A. There are 6 known binaries from high- 
resolution imaging in the sample, we count individual 
components if they fall into our selected spectral type 
range. None of the binaries show variability. In total, 
we find 19 late-L and 16 early-T dwarfs to match our 
criteria. Of those, 2 late-L and 5 early-T are variable at 
> 1% level with at le ast > 95% confidence. Following 
iBiirgasser et al.1 (j2003H we calculate the probability dis¬ 
tribution for the variability fraction for the late-L and 
early-T dwarfs using binomial statistics. We then find 
the combined probability for finding both a variable late- 
L and early-T dwarf to be 3.3liJ;3%. In reality, the prob¬ 
ability may be a bit higher because variability amplitudes 
are known to vary with time and may sometimes result 
in non-detections even if the object is variable at other 
times. This is clearly the case for Luhman 16A. 

Nevertheless, at a probability on the order of 10% or 
less, finding these large amplitudes in both objects may 
be more than chance. One possibility is the spin-axis: 
Luhman 16B is most likely seen near equator-on and any 
heterogeneities would therefore translate into the max¬ 
imum possible amplitude. The likelihood for B to be 
variable may therefore be higher than average. For Luh¬ 
man 16A the case is less clear, since its possible rotation 
period may indicate that its axis is somewhat misaligned 
from B. However, without a reliable period we cannot 
confidently conclude on the orientation of the spin axis. 

Other shared properties of the binary are the age and 
initial composition. Furthermore, the objects have rela¬ 


tively similar effective temperatures and for both of them 
non-equilibrium chemistry appears to be important (cf. 

While we currently do not have a large enough sam¬ 
ple to look for correlations between these factors and the 
occurrence of variability, it is possible that some of them 
contribute to a higher probability of variability. 

4.2. Comparison to other observations 

Because of an overlap in spectral coverage between 
the G141 and G102 observations, we can directly com¬ 
pare the variability amplitude obtained one year apart 
at least for a small wavelength range. For Luhman 16B, 
a variability amplitude of about 11% was found at I.l- 
1.15 /rm, slightly larger than in the J and H band peaks. 
We now find an amplitude of about 9% at 1.1-1.15 ^m, 
remaining equal to shorter wavelengths. Assuming that 
only the overall brightness level has changed, but not the 
spectral dependence of the variability, we can conclude 
that the variability amplitude is slightly larger in the Y 
band peak than in the J a nd H bands. This is con sistent 
with the results found in iBurgasser et al.l (|2014ll . who 
found a gradual decrease in the amplitude with wave¬ 
length across the near-IR, and t he unresolved simultane¬ 
ous photometry observations bv iBiller et al.l (j2013l l. On 
the other hand, sp atially resolved observations from the 
iBiller et’all (1201311 study obtained in a different night 
showed a large amplitude in H band but no variability 
in J band, which is very difficult to explain with current 
models. 

Interestingly, no significant variability (< 0.5%) was 
found for Luhman 16A over the full G141 wavelength 
range, including the 1.1-1.15 ^m region. In the new 
G102 observations however, the variability amplitude is 
about 4% in that wavelength region and at shorter wave¬ 
lengths. This suggests that Luhman 16A goes through 
periods of quiescence and periods of variability. Extrapo¬ 
lating the measured variability amplitude to longer wave¬ 
lengths using the spectral dependence of th e variable L 
dwarfs (2M1821, 2M1507 [Yang et ^1201511 would sug¬ 
gest a variability amplitude of about 2-4 % in the J and 
H band. 

The fact that Luhman 16A is variable at a level of a few 
percent with a likely rotation period that is very similar 
to that of Luhman 16B, and the fact that its amplitude 
appears to change with time, may perhaps account for 
some of the long-term changes in t he light curve s prev i- 
ously attributed to Luhman 16B bv iGillon et al.l (j2013H . 
Furthermore, the interpretation of the light curves is not 
as straight-forward as if one of the two objects was non¬ 
variable. For example, it cannot simply be assumed that 
the variability amplitude of Luhman 16B was twice that 
measured in the unresolved observations, as some peaks 
and valleys in the light curves of Luhman 16B may ei¬ 
ther be reduced or enhanced by the variability of Luh¬ 
man 16A. It should also serve as a caution to observa¬ 
tional programs that aim to measure the variability of 
Luhman 16B u sing Luhman 16A as n on-variable refer¬ 
ence, as done in iBurgasser et al.1 (j2014H . Although it may 
have been a valid assumption at that particular point in 
time, it is not possible to validate this. Nonetheless, 
with the amplitude of Luhman 16B being at least twice 
as large as for Luhman 16A, it will remain the dominant 
contributor to the combined variability. 

Our observations suggest the persistence or re- 
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emergence of a dark spot on Luhman 16B over the time 
scale of one year. It may be the same prominent dark 
spot that is also present in the D oppler imaging obser¬ 
vations by iCrossfield et al.l (|2014[ 1 in the K band. 

Spatially resolved time series of Luhman 16A and B 
at wavelengths < 1 u m have so far only been obtained 
bv iBiller et al.l (1201311 in the r’, i’ and z’ filters. The z’ 
filter covers approximat ely the same wavele ngth range as 
our G102 observations. IBiller et al.l (|2013ll found the z’ 
band light curves of Luhman 16B to be in phase with 
the J and H band light curves, this is consistent with our 
observations, since observations with the G141 did not 
show phase shifts between the 1.1-1.15 ^m region and the 
J and H bands. On the other hand, the i’ band light curve 
was anti-correlated with the z-band light curves. Our 
observations only reach the edge of the i’ band, but the 
steeply decreasing amplitude below 0.85 /xm may point 
to a phase reversal in the deep potassium feature that 
dominates the i’ band. On the other hand, the tentative 
detection of i’ and z’ band variability in Luhman 16A had 
these light curves in phase, with amplitudes of 3% in the 
z’ band and 2% in the i’ band. However, we observe a 
similar drop in the amplitude below 0.85 /xm, from about 
4 % to below 2%. Considering the very similar variability 
behavior across the 0.8-1.1 /xm region for Luhman 16A 
and B, the fact that only Luhman 16B shows a phase 
reversal in the i’ band is somewhat surpri sing. 

A recent survey bv iHeinze et al.l (j2015fl has indicated 
that some T dwarfs may be much more highly variable in 
the optical than the near-infrared. They found relatively 
high amplitude optical variability in two T dwarfs for 
which lower-amplitude variability had previously been 
detected in the near-IR. However, because observations 
were not taken simultaneously, it is also possible that 
the amplitudes have simply varied in time. Our observa¬ 
tions suggest this to be the case for Luhman 1 6A. The 
simultaneous monitoring by IBiller et’aTI (|2013ll already 
suggested that Luhman 16B had similar optical and near- 
IR amplitudes, although an unexplained low amplitude 
was found in one night only in the J band. Thanks 
to the overlapping wavelength region in the G141 and 
G102 grism, our observations confirm that the ampli¬ 
tude between 0.8 and 1.6 /xm remains fairly constant, 
with the exception of a decrease in the 1.4 /xm water 
ban d. Simultane o us op tical and near-IR observations of 
the IHeinze et al.l (|2015[ 1 sample would therefore be cru¬ 
cial to understand whether some objects indeed show 
strong differences between these wavelengths. If so, it 
might point to a different origin of the variability than 
the thickness variations of the silic ate clouds we propos e 
in this paper (cf. Sect. 14.31) and in lBuenzli et al.l (1201511 . 
That might be the case in particular if the variability 
of these objects originates predominantly from the dark 
0.7-0.8 /xm region dominated by the strong potassium 
opacity, especially if little variability is present at 0.85- 
0.9 /xm where the flux is much higher. 

4.3. Comparison to atmospheric models 

We use the same models as in iBuenzli et al.l (120151 ) 
to determine whether they are able to also explain 
the observations at 0.8- 1 .15 /x m. They are based on 
the lAckerman fc MarlevI (1200111 cloud models and have 
been updated with new opacities, among them FeH 
(|Saumon et al.l[2012|l . The models and parameter range 


covered are described in IBuenzli et al.l (1201511 . In brief, 
two models are linearly combined. These can have differ¬ 
ent effective temperatures Tes.i and Tf,s ,2 and different 
cloud parameters /sed,i and fsed, 2 - A smaller fsed value 
corresponds to a vertically and optically thicker cloud. 
Gravity log g and the disequilibrium parameter K^z have 
to be equal for the two models. The weights in the lin¬ 
ear combination are set using the covering fraction ci of 
the warmer model on the hemisphere where the emitted 
flux is maximal, and the change in covering fraction Aci. 
Then, the combined flux Fjnax = ciFi -|- (1 — ci)F 2 and 
Tjnin = (Cl - AcijAi -I- (1 - (Ci - Aci))F2. 

For Luhman 16A, the non-variability in the G141 ob¬ 
servations and the near-IR spectrum had suggested a ho¬ 
mogeneous cloud layer with intermediate cloud thickness 
fsed = 2 and Te/f = 1200 K. Because Luhman 16A is 
now found to be clearly variable, we re-fit the models 
using the full parameter space and allowing for two dif¬ 
ferent surfaces as for Luhman 16B. We find two potential 
solutions that are a good ht to the average spectrum and 
show a relatively flat variability amplitude between 0.9 
and 1.1 /xm as well as a sharp decrease below 0.85 /xm. 
These are shown in Figure [3 Both best models for Luh¬ 
man 16A have log// = 4.5 and K^z = 10"*’ cm^s“*^. The 
best-ht model remains the model where one hemisphere 
is fully homogeneous with fsed = 2 and Te/f = 1200 K, 
while on the other hemisphere 6% coverage fraction with 
cooler, thicker clouds {fsed = 1 and Teff = 1100 K) is 
introduced. An alternative possibility is a model where 
the cool, thick clouds {fsed = 1 and Teff = 1100 K) 
dominate with 82% coverage fraction, with the remain¬ 
ing 18% being warmer, thinner clouds {fsed = 3 and 
Teff = 1300 K). In this case, the change in thick cloud 
coverage fraction from one hemisphere to the other is 
only 1.5% to obtain a variability amplitude of about 4%. 
This model also results in a larger radius, 1.1 i?j, as com¬ 
pared to 0.95 Rj. The larger r adius is more in line wit h 
the expected evolution radius (|Saumon fc Marle3l2008ll . 
However, for this model, an amplitude decrease of about 
0.5% is predicted for the GrH band, and a decrease of 
1% below 0.9 /xm, which is not seen in the observations. 
We note that even thinner clouds or fully non-cloudy re¬ 
gions would result in even stronger amplitude differences 
inside and outside of absorption bands, while more thick 
or cooler clouds would make the spectrum too red. 

Overall, the best-ht models indicate that the near-IR 
variability of Luhman 16A, when present, would show a 
similar wavelength dependence as that of Luhman 16B, 
with a decrease of the variability amplitude in the water 
band and a slight gradient across the J and H band. 
However, the two Luhman 16A models could be clearly 
distinguished if variability could be measured also at 1.1- 
1.7 /xm and/or in the K band. The variability amplitude 
is predicted to be signiHcantly lower in the 1.4 /xm water 
band and K band for the more heterogeneous model with 
a larger fraction of cool, thick clouds. 

For Luhman 16B, the spectral dependence of the vari¬ 
ability in particular across the 1.4 /xm water feature had 
suggested a mixture of thin, warm {fsed = 3 and Teff = 
1300 K) and thick, cool {fsed = 1 and Tgff = 1000 K or 
1100 K) clouds, where the coverage fraction was some¬ 
what degenerate with other parameters such as the ra¬ 
dius and surface gravity. To fit the new observations 
and the overall spectrum, we Hnd that we require the 
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Figure 5. Comparison of observations and models for Luhman 16A. The two best 2-component models (blue and red) are shown with 
corresponding model parameters given in the Figure. Top: The maximum (dark blue and red) and minimum spectrum (light blue and 
orange), i.e. as given by the different covering fractions on the two hemispheres, for each model compared to observations. For clarity, 
the observations and models for one of the model cases has a vertical offset of 2 x erg/s/cm/A. The observations are shown 

in different shades of gray: The new G102 obs ervations from 0.8 t o 1.15 fim are shown in black (maximum spectrum) and dark gray 

1 minimum spectrum). The G141 (1.1-1.65 um. IBuenzli et ’^1201^) observations are shown in a medium gray and the FIRE spectrum 
Burgasser et al.ll2Q131 in light gray. The ground-based FIRE observations are strongly impacted by telluric lines in the water absorption 
bands. Bottom: The relative amplitude as a function of wavelength, i.e., the difference of the maximum to the minimum divided by the 
mean, for the observations (black) compared to the two models (red and blue), and model predictions for longer wavelengths. 


same overall model parameters except for different rela- 
tive coverage fractions (Fig. IH]). The coverage fraction 
Cl (fraction of the warmer, thinner cloud on the hemi¬ 
sphere with maximum flux) is reduced by 3-7% to 53 or 
78% coverage with the thinner clouds, and the difference 
in coverage fraction Aci between the two hemispheres 
is reduced by 0.7-2% to 5.3 or 8%, depending on the 
model. Both models result in a very flat spectral depen¬ 
dence of the amplitude, including across the FeH feature, 
very similar to what was observed. The similarity of the 
model results for the two different wavelengths and times 
suggest that it is primarily the spatial distribution of the 
thin and thick clouds that changes in time, and not the 
intrinsic cloud properties of the two components. 

The steep decrease in the amplitude below 0.84 /rm is 
better represented by the higher surface gravity model, 
but the model fits cannot conclusively determine the sur¬ 
face gravity of the object. The overall fit of the mean 
spectrum between 0.8 and 1.15 /rm is not optimal. In 


particular, the spectral slope is underestimated by the 
models, and a few absorption features are predicted to 
be deeper than measured, especially the CrH feature at 
1.0 /rm. As is the case for Luhman 16A, introducing 
fully non-cloudy regions would result in much more vari¬ 
able amplitudes as a function of wavelength. That case 
had already been shown unfeasible to explain the spec- 
troscopi c variability measur ed in the G141 grism at 1.1- 
1.7 p,m (|Buenzli et al.l[2M^ . 

The modeling results suggest that the cloud structure 
of Luhman 16A and B may be very similar in terms of 
cloud thickness and temperatures, as both can be fit by 
model combinations of fsed = 1 and 3 and Teg = 1100 
and 1300 K. The main difference in that case is that 
the fractional coverage of thicker, cooler clouds domi¬ 
nates on Luhman 16A, while that of the thinner, warmer 
clouds dominates on Luhman 16B. It would suggest that 
at least in the early L/T transition from late-L to early- 
T dwarf, the spectral changes could be explained by in- 
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Figure 6. Comparison of observations and models for Luhman 16B. Two best 2-component models (blue and red) are shown with 
corresponding model parameters given in the Figure. Top: The maximum (dark blue and red) and minimum spectrum (light blue and light 
red), i.e. as given by the different covering fractions on the two hemispheres, for each model compared to observations. For clarity, the 
observations and models for one of the cases has a vertical offset of 2 x 10“^® erg/s/cm/A. The observations are shown in different shades of 
gray: The new G10 2 observations from 0.8 to 1.15 /i.m are shown in black (maximum spectrum) and dark g ray (minimum spectru m). The 
G141 (1.1-1.65 um. IBuenzli et ^12015(1 observations are shown in a medium gray and the FIRE spectrum HBurgasser et al1{2013l) in light 
gray. The ground-based FIRE observations are strongly impacted by telluric lines in the water absorption bands. Bottom: The relative 
amplitude as a function of wavelength, i.e., the difference of the maximum to the minimum divided by the mean, for the observations 
(black) compared to the two models (red and blue), and model predictions for longer wavelengths. 


creasing the fractional coverage of thinner clouds. This 
is similar to the incre asing fractional cove rage of cloud 
holes as suggested in iMarlev et al.l (|2010H . except that 
having remaining thin clouds instead of deep holes ex¬ 
plains the spectrally flat variability ampliudes and keeps 
the FeH strength equal. The periods of non-variability 
could be explained by times when the coverage fraction 
of the two cloud types is more homogeneously distributed 
along both hemispheres, as differences of the order of 1% 
are needed to introduce an easily measurable variability 
amplitude. 

4.4. FeH as a probe of cloud patchiness? 

The FeH molecule contains the cloud-forming element 
Fe, and the condensation of Fe into a cloud dramatically 
decreases the abundance of FeH in the gas phase. The ob¬ 
served strengthening of the feature from late-L to mid-T 
spectral type has been suggested as evi dence that cloud 
holes appear at the L/T transition (|Burgasser et al.l 


I2002t iLodders fc Feelevll^OOGHMarlev et al.ll2010ll . This 
line of reasoning assumes that holes would allow flux 
from well below cloud base, where Fe is still present as 
gaseous FeH, to emerge. Spectra of alternatively clear 
and cloudy patches would then presumably show varia- 
tions in the FeH abs orption band depth. Alternatively, 
ICushing et al.l (|2008ll showed that a sequence of mod¬ 
els with increasing fsed value could also reproduce the 
strengthening of FeH. Therefore, even patchiness with 
different cloud thickness, i.e. different sedimentation ef¬ 
ficiency, could also be expected to result in variability in 
the FeH band strength. However we do not observe any 
differences in the variability amplitude within the FeH 
feature and outside, both in the late-L Luhman 16A and 
the early-T Luhman 16B. The spectral modeling of the 
variability for both sources (Sect. 14.31) is incompatible 
with cloud holes, but does suggest patchiness with dif¬ 
ferent cloud thickness. What are the implications of this 
result on the model of strengthening of the FeH band 
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across the L/T transition based on cloud patchiness? 

There are a few complications that can affect the FeH 
feature strength that we must consider. Even above 
the cloud base some FeH re mains in the gas in equi- 
librium with the condensate (|Visscher et al.l I20I0H , al¬ 
though FeH is mostly found below the cloud deck. Fur¬ 
thermore other opacity sources, notably gaseous K and 
CrH, are also important in this spectral region. So the 
visibility of the FeH band is not entirely modulated by 
the cloud but rather by the interplay of several factors. 
This complexity is evidenced by the fact that both Luh- 
man 16 A and B h ave the same FeH feature strength 
(|Fahertv et al.l[2Mll . despite Luhman 16B being being 
3 sub-types later and with a thinner cloud deck than 
Luhman 16 A. Although the scatter in observed FeH 
strengths for late-L and early-T dwarfs was too large 
to suggest d efinite strengthening a lready at those spec¬ 
tral types in iBurgasser et all (1200211 . and the sample was 
small, the co-eval nature of the binary helps to confirm 
that the bluer color and larger variability for an early-T 
compared to a late-L at the same age and initial compo¬ 
sition do not necessarily coincide with the strengthening 
of FeH or different variability in the FeH feature. 

To help illustrate the depth from which flux emerges 
from the atmosphere as a function of wavelength for dif¬ 
ferent scenarios, in Figure [7] we present brightness tem¬ 
perature spectra for the two cloud components of one of 
the Luhman 16B models shown in Figure[6l We also plot 
two cloudless models that were not able to fit the observa¬ 
tions, but illustrate how the reached depth changes when 
the cloud opacity is removed. Here brightness tempera¬ 
ture serves as a proxy of atmospheric depth from which 
photons at a given wavelength escape. In reality the con¬ 
tribution function for each wavelength spans a range of 
pressures around the level in the atmosphere where the 
local temperature is equal to the brightness temperature. 
The approximate locations of the model cloud base and 
cloud tops are shown. Although cloud top is arbitrarily 
denoted at r ~ 1, additional cloud opacity is found above 
these lines. 

Figure [7] makes clear that flux from the model which 
comprises most of the composite flux shown in Figure [6l 
fsed = 3, Teff = 1300 K, emerges from above the nominal 
cloud top and well above the Fe cloud base. For the thick 
cloud component, fsed = 1, Tes = HOOK, flux emerges 
from just below the Fe cloud top, but still far above cloud 
base. Even for a comparable cloudless model (here with 
Tes = 1200 K), which highlights the contribution of the 
gas opacity in the atmosphere, the high molecular opac¬ 
ity at 1 ^m prevents probing to the typical cloud base 
levels. Only in a much hotter (1600 K) cloudless model 
would the gas opacity allow emergent flux at 1 ^m to 
reflect atmospheric conditions below cloud base. This 
figure leads us to conclude that for objects as cool as 
Luhman 16 A and B even the complete absence of Fe 
cloud opacity would not result in an FeH absorption fea¬ 
ture from the deep atmosphere to be fully resurgent. 

This conclusion, which is contrary to earlier sugges¬ 
tions for using FeH as a tracer of cloud holes or patch¬ 
iness, stems mostly from a better understanding of the 
alkali opacities in the same wavelength region and a bet¬ 
ter atmospheric models incor porating cloud opacity than 
were available at the time of iBurgasser et akl (l2002f) and 
iCushing et al.l (I2008D . Certainly for higher effective tern- 



Figure 7. Brightness temperature as a function of wavelength 
for different atmospheric models for Luhman 16B, as labeled along 
the right hand side. Numbers refer to model effective temperature 
while ‘nc’ denotes cloudless models; fl and f3 refer to models with 
cloud sedimentation efficiency /sed=l 3, respectively). All 
models assume logg = 5. The two cloudy models are the same as 
those that comprise the cloudy case shown in the bottom right of 
Figure 0 For the two cloudy models the approximate bases and 
tops (r 1) of each model’s iron cloud are shown by the horizontal 
solid and dashed lines, respectively. The top of each model’s silicate 
cloud is also shown as a dashed line, for these models the silicate 
cloud base is at a similar temperature to the iron cloud base. The 
FeH band is also indicated by the vertical pair of lines. 

peratures, particularly early- to mid-L variable dwarfs, 
the band strength of the 1 fim FeH band may well be 
variable. 

For mid-T dwarfs, where clouds appear to have mostly 
cleared out from the photosphere, the situation may be 
different than at the onset of the L/T transition. There, 
the strengthening of F eH is more robustly observed by 
iBurgasser et al.l (120021 1. even though we are most likely 
not able to see into potential holes in the Fe clouds. Since 
the silicate cloud top is predicted by our cloud model to 
overlie the iron cloud top, the silicate cloud may, counter¬ 
intuitively, be of greater importance than the Fe cloud to 
the FeH band depth. As Figure [7] makes clear, the tops 
of the combined silicate and iron clouds greatly influence 
the depth from which flux near 1 /im emerges from the 
atmosphere. As the cloud opacity in the photosphere 
drops significantly towards the end of the L/T transition 
(/sed ^ 4), more of the remaining FeH gas may become 
visible. Such a gradual sinking of th e cloud was initially 
ruled out bv IBurgasser et aH ^2002( 1 under the assump¬ 
tion that the FeH was depleted above the cloud. 

Alternatively, variations in the cloud top altitude along 
with the local temperature at the cloud top could result 
in variations in the continuum flux and might also al¬ 
low some variation in the FeH band depth if clouds are 
sufficiently thin. No conclusive variability has been ob¬ 
served for these mid-T dwarfs, although there were small 
trends observed for the T4.5 2MASS J055 91914-1404488 
(|Buenzli et al.l IMTl: iRadigan et al.l [2^111 . which has a 
significant FeH feature. We note that a heterogeneous 
cloud distribution does not necessarily have to result in 
variability if the distribution is rotationally symmetric. 
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distributed in many small patches, or if the object is seen 
close to pole-on. 

At the onset of the L/T transition, a spatially hetero¬ 
geneous decrease in fsed between Luhman 16A and B ap¬ 
pears to be responsible for the spectral changes and the 
observed variability, but the remaining clouds are still too 
thick to affect the FeH strength. It may be reasonable to 
assume that a similar combination of both mechanisms 
happens toward the end of the transition as well. Only 
more sophisticated modeling accounting for the columns 
of gaseous FeH above the cloud base and variations in 
the cloud top altitudes could rigorously test these con¬ 
ceptual models for the re-emergence of the FeH feature 
at mid-T type. Such models must account for possible 
heterogeneities in the different cloud layers, perhaps with 
varying horizontal distributions. Such a detailed study is 
well beyond the scope of this paper but provides fruitful 
possibilities for future research. 

For Luhman 16A and B, the wavelength dependence 
of the variability amplitude, i.e. the different amplitude 
in the water band but otherwise overall relatively gray 
variability, strongly points toward cloud top variations in 
the silicate cloud layer. Models with cloud holes, i.e. a 
mixture of a cloudy and clear photosphere, clearly cannot 
explain the observed variability, while a mixture of thin¬ 
ner and thicker clouds, i.e. cloud top variations, can. In 
the 1.4 fim water band the flux emerges from well above 
the silicate cloud top. The reduced variability in that wa¬ 
ter band is therefore fully consistent with the thin/thick 
model. The model also predicts the same light curve 
phase at all wavelengths except in the deep K I line at 
0.77 /im, while anti-correlation in the water band with 
respect to the J band peak might occur for some combi¬ 
nations of cloudy/clear models. 

A ^ 7% flux variation at 1 ^m can be produced by a 
cloud top temperature difference of about 14 K, just by 
the change in the Planck function. Given the local ther¬ 
mal profile, this corresponds to only a ~ 0.1 bar variation 
in the cloud top pressure. Such a scenario is also broadly 
consistent with the observed differences in the variability 
amplitude in the very deep K I line, which is extremely 
sensitive to pressure. 

Finally, there is also the possibility that the observed 
variability could stem so lely from thermal variations 
l|Robinson fc Marlevll20l4l rather than cloud height vari¬ 
ations. However, in this case there is not an obvious 
mechanism to maintain such a temperature anomaly, and 
mixed cloud properties are needed to reproduce the aver¬ 
age spectrum of Luhman 16B. Therefore, our new obser¬ 
vations further support the theory of silicate cloud height 
variations as the dominant reason for the variability we 
observe between 0.8 and 1.7 p,m. 

5. CONCLUSIONS 


We have obtained the first observations of spectro¬ 
scopic variability of any brown dwarf in the 0.8-1.15 /im 
range, covering the FeH feature, by observing the L/T 
transition binary Luhman 16 with HST/WFC3. Both 
brown dwarfs display significant variability with likely 
periods of about 5 h, even Luhman 16A, which was previ¬ 
ously found non-variable in the near-IR and in a par tially 
overlapping wavelength region (iBuenzli et al.ll201^. an d 
tentatively variable at i and z’ bands ( Biller et al.ll2013ll . 
With both objects variable at least some of the time. 


the interpretation of spatially unresolved observations 
is more complicated than previously thought, and Luh¬ 
man 16A cannot necessarily be used as a non-variable 
reference in spatially resolved observations. Finding two 
high-amplitude variables in a binary is somewhat un¬ 
usual, as the independent probabilities for a late-L and 
early-T to be variable are only on the order of about 
10%. The Luhman 16 system may therefore have a com¬ 
mon property that favors variability, although it is yet 
unclear what that common factor is. 

We show for the first time that the wavelength depen¬ 
dence of the variability amplitude from 0.8 to 1.15 /tm 
is largely flat. This is the case for both objects. It can 
be explained by cloud height variations in the silicate 
clouds, while we can exclude holes. A similar conclusion 
was already rea ched for two othe r early T dwarfs from 
color variations (lApai et al.ll2013ll and for Luhman 16B 
from the 1.1-1.65 /im observations (|Buenzli et al.llMTsIl . 
but here we can for the first time fit the 0.8-1.65 /im 
wavelength range covering many atomic and molecular 
features and obtain a consistent result. It suggests that 
the spectral changes across the L/T transition, at least 
from late-L to early-T, do not stem from the opening of 
deep holes, but from cloud thinning that can be spatially 
heterogeneous. 

We have also re-evaluated the role of FeH as a t r acer o f 
cloud holes, as first suggested bv iBurgasser et all (1200211 . 
or as a tracer of ch anging fsed value as suggested by 
ICushing et al.l (|2008ll . The models using the most recent 
opacities suggest that even models with a clear photo¬ 
sphere at 1200 K reach only the top of the Fe cloud 
at 1 /tm and observations therefore cannot constrain 
whether there are holes in the Fe cloud. However, this 
also suggests that the re-emergence of the FeH feature 
through the L/T transition does not stem from holes in 
the Fe cloud. The reason for the re-emergence of FeH, 
in particular for the mid-T dwarfs, is therefore still not 
fully solved. Vertical mixi ng was already excluded by 
iBurgasser et al.l (1200211 and lLodders fc FeglevI (1200611 . It 
is still possible that a decrease in cloud opacity or varia¬ 
tions in the cloud height and temperature in the silicate 
and/or iron cloud might affect the depth of the FeH fea¬ 
ture for mid-Ts due to more of the FeH gas above the 
Fe cloud becoming visible. Detailed model calculations 
are required that study the effect of the different cloud 
layers and include a detailed FeH abundance profile be¬ 
tween the clouds. However, for our late-L and early-T 
dwarf, where the cloud opacity is still very significant, 
the variability is dominated by cloud-height variations 
in the silicate cloud, and these are not sufficiently large 
to also affect the FeH strength. The fact that Luhman 
16 A and B have ap proximately the same FeH strength 
(|Fahertv et al.l[2M^ . do not show any variability in FeH, 
but certainly have patchy clouds of different thickness, 
clearly indicates that the changing FeH strength does 
not correlate with other cloud evolution tracers at the 
early L/T transition. 

Nevertheless, our observations have provided very 
valuable input into the structure of the silicate clouds 
at the transition from L to T-type brown dwarf. With 
their unique brightness, Luhman 16 A and B continue 
to be the prime target for detailed atmospheric stud¬ 
ies of L/T transition brown dwarfs. The fact that both 
objects show variability at least some of the time pro- 
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vides a window into the three-dimensional atmospheric 
structure for two co-eval objects. The James Webb Space 
Telescope (JWST) will be able to precisely measure their 
variability over a broader wavelength range and at higher 
spectral resolution. Although the projected separation of 
the binary is currently decreasing and it may soon not 
be possible to spatially resolve them, it is expected that 
the projected separation will be increasing again by the 
time JWST launches and will be large enough to spa- 
tially resolve the binary . A po tential orbit suggested by 
iSahlmann fc Lazorcnkd (120151 1 suggests closest approach 
in 2017 and a separation of > 0.5" by 2018, increasing 
rapidly thereafter. 

Furthermore, astrometric monitoring will yield dy¬ 
namical masses that can be used as independent 
inputs into atmosphe r ic an d evolutionary models. 
Sahlmann &: Lazorenkd ()2015r i already provide a first es¬ 
timate of the system mass ratio of q — 0.78±0.10, con¬ 
firming the somewhat lower mass of the further evolved 
B companion. However, the luminositi es of the two 
brown dwarfs appear to be nearly equal (iFahertv et al.l 
12014 iLodieu et al.l I2015fl . The system therefore likely 
has a very shallow mass-luminosity relation, similar to 
the L/T transition binary SDSS J105213.51-b442255.7 
(jPuDuv et al.ll2015ll . This lends further support to the 
necessity for evolutionary models to take into account 
cloud evolution at the L/T transiti o n, as done for the 
hybrid tracks of ISaumon fc MarlevI (120081) . Individual 
masses will be determined once about half of the bi¬ 
nary orbit is covered. The orbit is likely highly in- 
clined and a full orbital perio d may be about 45 years 
(jSahlmann fc Lazorenkd 120151) . Although it may take 
about 20 years, the Luhman 16AB system will certainly 
become an extremely valuable benchmark L/T transition 
binary. Together with the variability, which we have now 
shown to be present for both objects and across a broad 
wavelength range, it will allow for detailed testing of at¬ 
mospheric and evolutionary models that take into ac¬ 
count the heterogeneous 3D atmospheric structure and 
evolution. 
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